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Abstract 

We present a model in which, due to the quantum nature of the signals 
controlling the implementation time of successive unitary computational 
steps, physical irreversibility appears in the execution of a logically re- 
versible computation. 

Introduction 

Quantum computation is mostly seen as a static process: we are given an al- 
gorithm A that must be decomposed in elementary computational steps which 
have to be unitary operators acting on the state space of the input /output reg- 
ister; the computing process consists then in the sequential application of the 
computational primitives to the register. This approach is the same adopted in 
classic computer science and has its genuine motivation in the functioning of a 
Turing machine in its simplest version: at each time step the tape is read and 
an operation is performed. In our desktops the time steps correspond to the 
ticks of the internal clock of the processor which "beats the bongos" at a given 
frequency. 

As soon as we leave the computational paradigm to come down to the imple- 
mentation of quantum algorithms, it turns out that the computational process 
is the time- evolution of a physical system prepared and acted upon by external 
agents in such a way that the initial state | i?(l) ) of the register evolves, at a 



certain time t, into the final state A\ R(l) ) obtained as 

A|JZ(1)) = U s ^U a _ 2 ...U 2 U x \R{l)) = 

= exp(-ii? s _iT s _i) exp(-iff s _2T s _ 2 ) • • ■ 
. . .exp(-ii/ 2 T 2 ) exp(-iifin)| R(l) ) 

Hi, i = 1, . . . , s — 1, being suitable Hamiltonian operators. This kind of external 
control (switching on and off the Hamiltonians at given times) is a physically 
very hard task to perform. The alternative definition of a single time indepen- 
dent Hamiltonian K under which the register can autonomously evolve from 
|-R(1)) toil|i?(l)) = exp(— iKt)\ R(l) ) turns out to be an equivalently difficult 
problem (see, for example, HDJ). 

Feynman's proposal for a quantum computer j^j represents a way out from this: 
by coupling the input/output register to additional degrees of freedom acting 
as a clock, it becomes immediate to define a single time independent Hamilto- 
nian determining the desired evolution. Furthermore, being able to implement 
the Toffoli gate, the computing model is able to compute all the function com- 
putable by a deterministic reversible Turing machine. 

Some of the features of the model have already been investigated in previous 
work In particular, in £j] we have shown that the use of conditional 

jumps makes Feynman's computer an implementation of the continuous time 
quantum walk computational paradigm 

In this work we analyze the coupling between the input /output register and the 
clocking mechanism and the appearance of physical irreversibility in the context 
of logically reversible computation. 

Feynman's cursor model 

Let p(l) = | R(l) )( R(l) |, with | R(l) ) 6 T~lregister, be the initial state of the 
input/output register. Let U\, U 2 , ■ ■ ■ , U s -\ be the unitary operators represent- 
ing the successive "primitive" steps of the computation to be performed. 
Set 

| R(x) ) = U a -x | R(x - 1) >, 1< x < s. (1) 

Call/Cr the subspace of Tt re gister spanned by the vectors \ R(1)}, \ R(2)), . . . , \R(s)) 
Set 

d = dimQC r ). (2) 

Following the approach of (SJ, we model the clocking mechanism, which se- 
quentially applies the transformations Ui, U 2 , ■ ■ ■ , U s -\ to the register, with a 
quantum mechanical system, the cursor. 

We call Ti-cursor the s— dimensional state space of this system and refer it to a 
selected orthonormal basis | C{1) ), | C(2) ),..., | C(s) ). 

We call position of the cursor the observable Q acting on the vectors of this 
basis as 

Q\C{x)) =x \C(x)}, l<x<s. (3) 
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We suppose that the state of the overall system, the machine, evolves in the 
HUbert space r H ma chine = ^register <8> Hcursor under the action of a Hamiltonian 
of the form 

A s_1 

H = --Y J U x ®\C{x + l)){ C{x) | + U" 1 <g> | C(x) )( C(x + 1) |. (4) 

cc=l 

We suppose, furthermore, that the state of the machine is, at time t — 0, 
represented by the vector 

|M(0)) = |i2(l))®|C(l)>. (5) 

It is well known [7] that for every time t it is then 

s 

I M(t) ) = J2 *5 s) | fl(x) ) ® I C(a:) ) (6) 

where 

2 s 

c(t, x; s) = } exp [iXt cos(i9(fc; s))l sin(t?(fc; s)) sin(a; s)) (7) 

s + 1 ti 

and 

tf(fc; S ) = A?l (8) 

Call 

p m (t) = \M(t))(M(t)\ (9) 
the density matrix of the machine at time t. 

By taking the partial trace Trn eurBor {pm{t)) with respect to the cursor degrees 
of freedom, we get the density matrix p r {t) of the register: 

s 

p r (t) = J2Ht^:S)\ 2 p(x), (10) 
x=l 

where, for 1 < x < s, 

p(x) = U x -i . . . U lP (l)U^ . . . U~\. (11) 

In order to trace out the register degrees of freedom, it is expedient to refer IC r , 
at each time t, to the orthonormal basis | b\{t) ), | b^it) ), • • • , | b^{t) ) formed by 
the eigenvectors of p r {t): 

Pr{t)\b i {t)) = \ j {t)\b i (t)), l<j<d. (12) 

A simple computation shows, then, that the density matrix of the cursor is given 

by 

Po(t) = J2 X j( t )\d j (t))(dj(t)\ (13) 
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where the sum extends to the values of 1 < j < d such that Xj(t) > and 
where, for such values of j, | dj(t) ) is given by 

| dj(t) > = -j== J2 <t, x; s)( bj(t) | R(x) ) | C{x) }. (14) 

As it is easy to check that ( dj(t) \ dk(t) ) = Sj.k, the above representation of 
p c {t) shows that it has eigenvalues Xj(t) and, therefore, von Neumann entropy 

S(Pc(t)) = -'£ i \ j (t)]n\ j {t) = S(p r (t)). (15) 

3 

The equality, in the particular case considered here, in which the initial state 
of the machine is the pure state | M(0) ) = | -R(l) ) ® | C(l) ), between the 
entropy S(p c (t)) of the cursor and the entropy S(p r (t)) of the register (a well 
known property of a bipartite system in a pure state) is easily understood, in 
physical terms, by the insertion, in the expression for | M(t) ), of the partition 
I r = I bj(t) ) ( bj(t) | of the identity in JC r : 

d s 

3 = 1 x=l 
d s 

= E I b ^ ) ® E <^ s )< h i (*) I R W ) I C{x) ) = 

d / — 

= Ev A ^)i^(*)>®i^(*)>- (16) 

The above expression (the Schmidt decomposition of the state | M(t) )) shows 
that if upon a measurement at time t, the register is found in state | bj(t)), then 
the cursor collapses into the state | dj (t) ) , and vice versa. 



An explicitly solvable example 

We focus our attention, in what follows, on the simplest non trivial case, in 

which dim(Ti reg i S ter) — 2. We consider, namely, the simple case in which the 

register is a two level system or, equivalently, a spin 1/2 system. 

We indicate by a — (a\, oi, o^) the three components of such a spin in an 

assigned reference frame and by e^e^e^ the versors of the three coordinate 

axes. 

In the basis | o 3 = ±1 ), the density operator p r (t) will be represented by a 
matrix of the form 

1/ l + s 3 (t) s 1 (t)-is 2 (t) 
Pr[, ~2\s 1 (t)+is 2 (t) l-« 3 (t) 
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where 

s j (t)=Tr(p r (t)-a j ), j = 1,2,3. (18) 

Equivalently stated, the Block representative of the state p r (t) is given by the 
three-dimensional real vector 

s(t) = si(t) ■ e x + s 2 (t) ■ e 2 + s 3 (i) • §3 = 

S 

= 5>(i,x; S )| 2 Tr(p(x)-a). (19) 

x=l 

We examine, here, the behaviour of s(t) in the simple example defined by the 
following additional conditions: 

i. The initial state of the register is: 



I R(l) > - COS ^- j I (7 3 = +1 ) + sin I - 1 I (73 = -1 ), (20) 

namely the eigenstate belonging to the eigenvalue +1 of n(l) - a, with 

71(1) = e x sin# + e 3 cos 6*. (21) 

ii. Each of the unitary transformations U x is a rotation of a fixed angle a 
around the axis e 2 

<7i = [/2 = ... = t/ s -i = e~^' T2 . (22) 

We wish to remark that the above example captures the geometric aspects not 
only of such simple computational tasks as NOT or V NOT (viewed as rotations 
of an angle it or it/ 2 respectively, decomposed into smaller steps of amplitude 
a) but also of Grover's quantum search |S]- 

If the positive integer fi is the length of the marked binary word to be retrieved, 
set 

x(/i) = arcsin(2~^) (23) 

and 

6 = 7T-2 X (p). (24) 

Then the state (12011 correctly describes the initial state | t ) of the quantum 
search as having a component 2 - ' 1 / 2 in the direction of the target state, here 
indicated by | co ) — \ (73 — +1 ), and a component Vl — 2~f* in the direction of 
the flat superposition, here indicated by | (73 = — 1 ), of the 2 M — 1 basis vectors 
orthogonal to the target state. In this notations, if 

a = -4 X (/i), (25) 

then the unitary transformation exp(— iaaij2) corresponds to the product B- A 
of the oracle step 

A = I r - 2 I w ) { uj I (26) 
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and the estimation step 

B = 2\ l)(l\- I r (27) 
where I r is the identity operator in Ti re gister ■ 

It is having in mind the connection with Grover's algorithm that, for the sake 
of definiteness, in the examples that follow we are going to consider the one- 
parameter family of models, parametrized by the positive integers /x, corre- 
sponding to the choice l|24[l and (|25(l of the parameters 9 and a and to the 
choice s = 2^ + 1 of the number of cursor sites, corresponding to the possibility 
of performing up to an exhaustive search. 
In the example defined by (|20|l and 122|) it is 

( R{x) | a\ R(x) ) = sin (9 + (x - l)a) e x + cos (9 + (x - l)a) e 3 (28) 

and, therefore, 

s 

s(t) = ^2 l c (*> ^ s )| 2 ( sin {0+(x- l)a) e t + cos (9 + (x - l)a) e 3 ) . (29) 

x=l 

Figure^presents a parametric plot of (si(i), S3(i)) under the above assumptions. 




Figure 1: A parametric plot of (si(t), ss(t)) for < t < s, A = 1. The choice 
/i = 7, x — arcsin(l/2 M / 2 ), s = 2 M + 1, a = —4%, 9 = tt — 2\ of the parameters 
is motivated by the connection with Grover's algorithm. 

It is convenient to describe the Bloch vector s(t) = si^)^ +S3(i)e 3 in polar 
coordinates as 

si(i) =r(t)smj(t), s 3 (t) = r(t) cos j(t). (30) 
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The eigenvalues of p r (t) can then be written, in this notation, as 

A 1 (t) = |(l+r(«)) ) A 2 (i) = i(l-r(<)) (31) 
and the corresponding eigenvectors as 

\h (+\ \ — ( cos (7(0/2) \ \h (t\ \ — ( ~ sin(7W/2) \ 

1 bl{ - l) > ~ { sin( 7 (t)/2) ) ' 1 h{t > > - { C0B( 7 (i)/2) J ■ (32) 

The two corresponding states of the cursor, in the sense of (|16ll . are then: 
\d 1 (t)) = -^L=J2c(t,x;s) cos((0+(x-l)a-7(t)/2)) | C(i) ); (33) 

\d 2 (t)) = -^=J2c(t,x;s) sin((0+Or-l)a- 7 (i)/2)) | C(i) ). (34) 

The von Neumann entropy 5 (p r (t)) is therefore 

S(pr(t)) = 2 ln 2^ 2 2 ' ' 

An example of its behaviour is shown in figure [3 It is to be stressed that, as 



50 100 150 200 250 

Figure 2: The von Neumann entropy of the register as a function of time, for 
the same model as in figure ^ for < t < 2s. The dashed part of the graph 
shows the effect of the cursor wave packet being reflected at the rightmost site 
s. 

Ai(i) > A 2 (i), at each time t the projector |&i(f))(&i(t) | is, among the projectors 
on the state space of the register, the one having in the state p r (t) the greatest 
probability of assuming, under measurement, the value +1. 
Stated otherwise, if the desired output of the computation is a given state | B\ ) 
(the eigenstate corresponding to the eigenvalue +1 of the component of a along 
an assigned direction in the e^eg plane) the optimal choice of the time t at 
which to read the output is such that | b\ (t) ) = \B\). 

What figures Q and [21 show is that for no choice of t > is the probability of 
finding the "target" output \Bi) equal to 1: it is bounded above by Ai(i). There 
is always, as shown in figure El a non vanishing probability (bounded below by 
A 2 (i)) of finding the register in the orthogonal, "undcsired", state | B% ). 
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Figure 3: The thick line is a graph, for < t < s of Tr{p r {t){I + cr 3 )/2), 
the probability of observing the "target" state | ui ) = \ o% = +1 ) = | B x ) 
in the example of Grover's algorithm. We have evidenced the instant r = 
0(2 /i/ ' 2 ) at which this probability reaches its maximum. The dashed line is 
a graph of Tr(p r (t)(I — cr 3 )/2), the probability of observing the "undesired" 
output | <t 3 = -1 ) = | B 2 ). 

Reading the register 

With reference, for definiteness, to the example of figure 3, call r the instant of 
time at which the probability Tr(p r (t) \u>)(u> \) reaches its first and absolute 
maximum. We recall that, in the above example, the target state | to ) is taken 
to be the "up" state | a% = +1 ) of the register. 

The whole point of the analysis of the previous section is that Ai (r) is strictly 
smaller than 1 . This amounts, in turn, to a deficit 1 — Ai(r) in the probability 
of finding the target state. This deficit is not, in itself, a strong limitation in a 
quantum search algorithm, because we can in principle identify the right target 
through a majority vote among a "gas" of a large number TV of machines. The 
trouble is that if we want to use the same machines once more, we need to purify 
the "gas" of registers from the fraction A 2 (r) of them which have collapsed into 
the wrong state: standard thermodynamic reasoning 9 shows then that this 
requires the removal from the gas, supposing a heat reservoir at temperature 
T is available, of an amount of heat of NkBTS(p r (T)), ks being Boltzmann's 
constant. 

We wish, in this section, to supplement the above considerations with an ex- 
plicit description of the post-measurement state of the machine, showing, in 
particular, the effect onto the clock of the act of reading the register. 

Suppose that at the optimally chosen instant r, at which it is j(t) — 0, while 
the machine is in the state | M(r) ), a measurement of the projector (I r + 03) /2 
is performed. 

If the measurement gives the result 1, then the state | M(r) ) collapses to 

1 s 

I M 1 (r) ) = I a 3 = +1} ® -== V c(r, x; s) cos((0 + (x- l)o)/2) | C(x) ) . (36) 
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If, instead, the measurement gives the result 0, then the state | M(t) ) collapses 
to 

I M 2 (r) ) = | a 3 = -1 ) ® - 7 J= £ c(r, x; a) sin((0 + (x - l)a)/2)\ C(x) ) . (37) 

V A 2(r) x=1 

Figures 4. a and 4.b show the probability distributions 

Pi(x,t) = |(c(r, a; ; S )cos((0+(x-l) a )/2))| 2 /Ai(r) (38) 
P 2 (ar,r) = \{c(T,x;s)sm{{e + (x-l)a)/2))\ 2 /A 2 (r) (39) 

of the observable Q (position of the cursor) in the states | Mi(r) ) and | M 2 (t) ), 
respectively. Figures 5. a and 5.b show the probability distributions of the ob- 





Figure 4: Figures (a) and (b) represent, for the same choice of parameters as in 
figure ^ respectively the probabilities P\{x, r) and P 2 (x,r) as functions of x. 

servable H (energy of the machine) in the states | Mi (r) ) and | M 2 (r) ) , respec- 
tively. The two energy distributions of figures 5 are easily derived from the fact 
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Figure 5: Figures (a) and (b) represent the probability distribution pi(Ek) and 
P2(Ek) of the energy H in the state | -Mi(t) ) and | M 2 (r) ) respectively. 



that the Hamiltonian H defined in (4) has the eigenvalues 
Ek — —A cos($(fc; s)), fc = 1, . . . , s; 



(40) 
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each doubly degenerate, an ortho-normal basis in the eigenspace belonging to 
the eigenvalue E k being given, for instance, by the two eigenvectors 

s 

I E k ; a 2 = ±1 ) = | (J2 = ±1 ) <8 X] v *( x ) e MTia{x - 1)/2)| C(x) ), (41) 

x=l 

where 

v k (x) = sj sin(x §{k; a)). (42) 

This leads to the explicit expressions 

Pj(Ek)= E \{M 3 {T)\E k] o 2 = i 1 )\\ 3 = 1,2. (43) 
n=±i 

Figures 4 and 5 show that a collection of identically prepared and indepen- 
dently evolving machines becomes in fact, under the operation of reading the 
register at time r, a mixture of two distinct "molecular" species, "1" (present 
in a concentration Ai(r) ), and "2" (present in a concentration A2(r)). In each 
of these two molecular species, the same "atomic" constituents have arranged 
themselves in a different geometrical shape (figures 4), with a different orienta- 
tion of the register spin (equations and 1)23), because of a different energy 
distribution (figures 5). 

Comparison with the distribution of H in the pre-measurement state | M(t) ), 
given in figure 6, shows that the presence of the impurities of type "2" is due 
to unusually intense exchanges of energy between the machine and the reading 
(measurement) apparatus. 
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Figure 6: In the state | M(t) ) the probability distribution of H is given by 
p{E k ) = {v k {x)f. 



Discussion and outlook 

We don't claim that (quantum) computation cannot be made reversible. We 
have simply pointed out one aspect in which a reversible machine is an ideal- 
ization; this idealization amounts to neglecting the back-reaction (figures [SJ of 
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the clocked subsystem onto the clock. 

In very concrete physical terms, from the dispersion relation E(p) = — cos(p) 
(see equation J1DJ) it is immediate to conclude that the distortion of the en- 
ergy spectrum depicted in figure [S]b corresponds to a decrease in the speed 
v = dE(p)/dp of the cursor. This recoil effect could of course be neglected if the 
clock consisted of, say, 10 23 atoms (we refer the interested reader to the huge 
literature on the limitations posed by quantum mechanics to the measurement 
and/or operational definition of space-time distances Od]), but might be of 
relevance for a machine scaled down to a molecular chain evolving under the 
sole effect of its initial condition not being an eigenstate of the Hamiltonian. 
Stated in equivalent terms, there is nothing wrong in the assumption of starting 
the computation in a pure state of the register: we have pointed out, however, 
that realizing this initial condition has a cost NkBTS(p r {r)) in terms of energy 
to be dumped into the environment in order to get rid of the entropy generated 
(for the simple fact that the two recoil patterns of figures EJa and 03b have both 
strictly positive probability) in the previous run of the machine. 
The toy model corresponding to the choice has allowed us to show in an 
explicit quantitative form the decoherence induced by the coupling with the 
timing apparatus, appearing through the build-up of entropy in the state of the 
register subsystem. 

Beyond the details of the model considered, it is to be stressed that such a 
build-up is a general consequence of the fact that the coefficients c(t, x\ s) ap- 
pearing in (|10|) are determined by the discretized version of the free Schrodinger 
equation 

* jI c (*' x; s) = -- (c(t, x - 1; s) + c(t, x + 1; s)) , (44) 

leading to the well known phenomenon of wave packet spreading (quadratic 
increase in time of the variance of Q). 

We are well aware that the explicit model discussed in the previous sections 
is, under many respects, far from being optimized from the point of view of 
minimizing the probability deficit 1 — \\{t). 

The initial condition (jSJ), for instance, corresponds to the classical intuition of 
initially placing the clock in a sharply defined position | C(l) ). It will be worth 
studying the effect of a better choice of initial conditions, with the probability 
amplitude of the cursor spread on an initial extended region; in much the same 
spirit, we recall the analysis leading in to the proposal of supplementing the 
active part of the cursor subsystem (the collection of sites x for which U x is 
different from the identity) with an extended inactive part, having the effect 
of temporarily lowering the entropy of the register subsystem. The problem of 
optimally investing, in the above two ways (to the left and to the right of the 
active region) , an assigned amount of space resources in order to minimize at a 
selected time r the deficit 1 — Ai(t) in the probability of finding the target state 
poses itself as a natural question in this context. 

A related problem is that of providing a stability analysis of motion under 
position dependent coupling constants X(x) in l|4^|) . 

Under two more respects, in studying a model of the class defined by (@J, we 
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have not fully exploited the potentialities of Feynman's approach: 

• in Feynman's full model the quantum walk performed by the cursor is by 
no means restricted to a linear graph: the use of conditional jumps allows, 
in fact, to explore much more interesting planar graphs (we refer to |3| for 
a quantitative study of these more general systems); 

• the notation adopted in Q does not give a full account of the original 

intuition (better described in terms of creation and annihilation operators 

r±(x) ) of a single particle in a quantum lattice gas jumping between 

nearest neighbor sites. 
Having written the Hamiltonian Q as 

A s_1 

H = --^2U X ® t + (x + 1)t_(x) + If- 1 <Z> t + (x)t_{x + 1), (45) 

x=l 

the idea emerges quite naturally of studying the evolution under (|45|l of many 
particles. 

Possible applications of this proposal of a "multi-hand quantum clock" , its use in 
the repeated application of a given transformation, the steering of the particles 
along different branches of the graph in order to act in parallel on distinct parts 
of the register, the solution of conflicts in the application of non commuting 
primitives to a same part of the register deserve, we think, further research. 
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